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A multi-orbital model for a CoCb-layer in NaccCoCb is derived. In this model the kinetic energy for 
the degenerate t2^-orbitals is given by indirect hopping over oxygen, leading naturally to the concept 
of four inter-penetrating Kagome lattices. Local Coulomb interaction couples the four lattices and 
an effective Hamiltonian for the interaction in the top band can be written in terms of fermionic 
operators with four different flavors. Focusing on charge and spin density instabilities, a big variety 
of possible metallic states with spontaneously broken symmetry are found. These states lead to 
different charge, orbital, spin and angular momentum ordering patterns. The strong superstructure 
formation at x = 0.5 is also discussed within this model. 



I. INTRODUCTION 

The layered Na^Cc-02 has been initially studied for its 
extraordinary thermoelectric properties and for its in- 
teresting dimensional crossover But recently wider 
attention has been triggered by the discovery of super- 
conductivity in hydrated Nao.35Co02 and the discovery 
of an insulating phase in Nao.5Co02^*2i£ Since then, 
various types of charge ordering phenomena in Na x Co02 
have been reported, 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 but 
also strong spin-fluctuations and spin density wave tran- 
sitions have been observed, 23 ! 24 ! 25 ! 26 ! 27 ! 28 ! 29 ! 30 ! 31 ! 32 ! 33 ! 34 

The material consists of Co02-layers where Co-ions are 
enclosed in edge-sharing O-octahedra. These layers alter- 
nate with the Na-ion layers with Na entering as Na 1+ and 
donating one electron each to the Co02-layer. 

The electronic properties are dominated by the 3d-t2 g 
electrons of the Co-ions which form a two-dimensional 
triangular lattice. However, the spatial arrangement of 
the Na 1+ -ions plays a crucial role too for the physics of 
this material. There are two basic positions for the Na- 
ions, one directly above or below a Co-site and another in 
a center position of a triangle spanned by the Co-lattice. 
The metallic properties are unusual and vary with the 
Na-concentration and arrangement. 

A brief overview of the present knowledge of the phase 
diagram of Na x Co02 leads to following still rough pic- 
ture. The most salient and robust feature, at first sight 
is the charge ordered phase for x = 0.5 separating the Na- 
poor from the Na-rich system. The Na-ions arrange in 
a certain pattern inducing an insulating magnetic phase 
in the Co02-layers below 50K. 33 On the Na-poor side 
(x < 0.5) the compound behaves like a paramagnetic 
metal. When it is intercalated with H2O superconductiv- 
ity appears between x « 0.25 and x « 0.35. In several re- 
spects more interesting is the Na-rich side where one finds 
a so-called Curie- Weiss metal. Here the magnetic suscep- 
tibility displays a pronounced Curie- Weiss-like behavior 
after subtracting an underlying temperature independent 
part: x = C/(T — 0) where ranges roughly between 
-50 and -200 K depending on x, and the Curie constant 



is consistent with a magnetic moment in the range of 1 
- 1.7 \i B- Note that deviations from the Curie- Weiss be- 
havior have been observed at low temperatures™ On the 
Na-rich side a transition at high temperature ~ 250 - 340 
K has been observed and interpreted as crystal structure 
or charge orderingi 5 ^! 2 ^. For Nao.75Co02 a magnetic 
transition occurs at 22 K and is most likely a commen- 
surate spin density wave or ferrimagnetic order which is 
rather soft towards magnetic polarizationjSLSLSSiSii In- 
terestingly this magnetic phase is metallic and has even 
a higher mobility than the non- magnetic phase. For Na- 
content x > 0.75 several magnetic transition at a simi- 
lar critical temperatures have been observed, but /iSR- 
data suggest rather an incommensurate spin density wave 

order 28 | 29 i 34 i 

The arrangement of the Na ions between the layers 
depends on the Na doping x and several superstruc- 
tures have been found^Siifi the clearest evidence for the 
superstructure formation is at x — 0.5 where the Na- 
ordering leads to a metal-insulator transition at low 
temperatures] 7 ! 8,12 But also away from x = 0.5, NMR ex- 
periments indicate the existence of non equivalent cobalt 
sites and phase separation™^ 

The complex interplay between Na- arrangement and 
the electronic properties poses a interesting problem. 
Various theoretical studies have mainly focused on single 
band models on the frustrated triangular lattice, in par- 
ticular in connection with the superconducting phase ig- 
noring Na-potentials. 36 ! 37 ! 38 ! 39 ! 40 ! 41 ! 42 There is also work 
done on multi-orbital models 43 ^ 4 ^ and density functional 
calculations have been performed^ 5 ^ 4 ^ 4 ^ 4 ^ 4 ^ 5 ^ 5 ^ 5 ^ Ac- 
cording to LDA calculations, the Fermi surface lies near 
the top of the 3d-^2^-bands. They form a large hole-like 
Fermi surface of predominantly a\ g character in agree- 
ment with ARPES experiments. 53,54,55 In addition the 
LDA calculations suggest, that smaller hole pockets with 
mixed a\ g and e' g character exist on the Y — K direction 
on the Na-poor side. 

At the T point the states with a\ g and e' g symmetry are 
clearly split, but on average over the entire Brillouin zone 
the mixing between a\ g and e' g is substantial. Koshibae 
and Maekawa argued that the splitting at the Y point 
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originates from the cobalt-oxygen hybridization rather 
than from a crystal field effect due to the distortion of 
the oxygen octahedra, because the crystal field effect in a 
simple ionic picture would lead to the opposite splitting 
of the dig and e' g states4£ There is also spectral evi- 
dence, that the low-energy excitations of Na^CoC^ have 
significant O 2p character. 56 Reproducing the LDA Fermi 
surface with a tight binding fit for the Co t2 g orbitals, it 
turns out that the direct overlap integral between the 
cobalt orbitals is much smaller than the indirect hopping 
integral over the oxygen 2p orbitals. 44 Therefore, it is rea- 
sonable to start with a three band tight binding model of 
degenerate t^g orbitals, where the only hopping processes 
are indirect hopping processes over intermediate oxygen 
orbitals. This approximation provides an interesting sys- 
tem of four independent and inter-penetrating Kagome 
lattices as it was already pointed out by Koshibae and 
Maekawa. 

Our study will be based on this model band structure 
which has a high symmetry. Within this model we ex- 
amine various forms of order that could be possible from 
onsite Coulomb interaction. The paper is organized as 
follows: In section [HJ the tight-binding model and the 
concepts of Kagome operators and pocket operators are 
introduced. In section IIIII an effective Hamiltonian for 
the local Coulomb interaction is derived and in section 
II VI this effective interaction is written in a diagonal form, 
by choosing an appropriate basis of SU(4) generators. In 
section the effects of small deviations from our simpli- 
fied tight-binding model are discussed. In section IVTl all 
possible charge and spin ordering patterns of our model 
and the corresponding phase transitions, are shortly de- 
scribed. In section I VIII the relevance of the above de- 
scribed collective degrees of freedom to Na x Co02 is dis- 
cussed by comparing the different coupling constants and 
by taking into account symmetry lowering effects. In sec- 
tion lVIlfl we apply our model to the Na-ordering observed 
at x — 0.5 and we summarize and conclude in section llXl 



II. TIGHT-BINDING MODEL 

We base our model on the assumption that the 3d- 
^-orbitals on the Co-ions are degenerate. Their elec- 
trons disperse only via 7r-hybridization with the inter- 
mediate oxygens occupying the surrounding octahedra 
(FigQ. As noticed by Koshibae and Maekawa the re- 
sulting electronic structure corresponds to a system of 
four decoupled equivalent electron systems of electrons 
hopping on a Kagome lattice 4 ^. The different sites, how- 
ever, are represented by different orbitals. Each of the 
three orbitals {d yzi d ZXl d xy } on a given site participate 
in one Kagome lattice, and the fourth Kagome lattice 
has a void on this site. The corresponding tight-binding 
model has the following form, 

fftb = E E £ k m ' C L^m>*, (1) 
krr mm' 




FIG. 1: Schematic figure of a C0O2 plane drawn with cubic 
unit cells. The edge-sharing of the oxygen octahedra around 
the Co ions is visualized. The edges of the cubes are oriented 
along the coordinate system (x,y,z). The triangular lattice 
of the cobalt is spanned by the vectors ai, a2, a3. (ai + a.2 = 
—as), a — |a*| is the lattice spacing. An oxygen 2p orbital and 
the cobalt ti g orbitals hybridizing with it by 7r-hybridization, 
are shown. 

where o^ ma . = e* k r cj m are the operators in 

momentum space of c\ mc7 which creates a ^-orbital 
(dy Zl d ZXl d xy ) with index m G {1,2,3} and spin a G {| 
, 1} on the cobalt-site r. N is the number of Co-sites in 
the lattice. 

(—[i 2tcos(ks) 2£cos(&2) \ 
2tcos(k 3 ) -fi 2£cos(/ci) , (2) 
2tcos(/c2) 2tcos(fci) —\i J 

with ki = k • a^, cf. Fig. [IJ The hopping parameter t = 
tp d / A > 0, where t p d is the hopping integral between 
the p y and the d xy or d yz orbital shown in Fig. A 
is the energy difference between the oxygen p and the 
Co-t2 g levels. The diagonalization of the matrix ek by a 
rotation matrix Ok G SO (3) 

Eorc m 'or' = %^ (3) 

mm' 

results in the three energy bands 



E^ = t + ty/l + 8 cos(fci) cos(A:2) cos(ks) — ji 

E\ = t- ty/T+ 8cos(&i) cos(/c 2 ) cos(/c 3 ) - fi (4) 

El = -2t - /1. 

These bands have the periodicity £^+b = > wnere 
the vectors Bj are defined by 

2tt 

ai • = — sin(0i - 0j) i,j G {1 ... 3} (5) 

with 0j = 27T//3. These three vectors Bj connect the T 
point with the three M points in the Brillouin zone (BZ), 
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and the vectors 2 Bj are primitive reciprocal lattice vec- 
tors. The bands of this tight-binding model have there- 
fore a higher periodicity than the bands of a more general 
model. This leads to the appearance of special symme- 
try lines (thin lines) and symmetry points (M' and K') in 
the Brillouin zone, shown in Fig. where the bands are 
plotted along the line r'-K'-M'-P. Within a reduced BZ, 
these bands correspond to the bands of a nearest neighbor 
tight-binding model on a Kagome lattice 43 . The density 
of states per spin and per reduced BZ is also shown in 
Fig. It has a logarithmic singularity at E = 2t and 
jumps from y/%/(2wb) to at E = 4t. 




FIG. 2: The original Brillouin zone (BZ) of the triangular 
lattice consists of four reduced BZs around of the F point (0) 
and the three M points (1,2,3). The symmetry points of the 
reduced BZs, M', K' and V' are symmetry points for the tight 
binding model in Eq. JIJ due to the higher periodicity of the 
bands. It is therefore sufficient to draw the bands along the 
lines r'-K'-M'-T'. The Fermi surface (FS) for x = 0.5 lies at 
E^ « 3.16£. The density of states per spin and per reduced 
Brillouin zone D is given in units of 1/t. It has a logarithmic 
singularity at E — 2t. 

The states, that are connected by the considered hop- 
ping processes form a Kagome lattice. Since in this way 
the Co02-plane consists of four independent and inter- 
penetrating Kagome lattices 43 , it is convenient to label 
the states belonging to the same Kagome lattice with an 
index I G {0,1,2,3}. This can be done with the vectors 

of Fig. [I] as 

a tirn ~ C R+a/+a rr) /m- (6) 

In this way, the operators aH m with fixed I create all 
the states off a Kagome lattice. In the following, these 
operators will be called Kagome operators. Their Fourier 
transform is given by 

n n - _2_ p iK.(R+a z +a™) U (? \ 
V ^ V R 

where the vectors k belong to the reduced BZ, labeled 
in Fig. [3 and r runs over the lattice spanned by the 
vectors 2a^. 



The BZ consists of four reduced BZs shown in Fig. 
An alternative labeling of the states is obtained therefore 
by defining the operators 

M _ p-iBj-a™ t (Q\ 
17 Km ° °K+Bjra> \ U J 

where the vectors Bj are defined in Eq. © and in addi- 
tion we set Bo = 0. As shown in Eq. (|A1|) of Appendix 
A, the transformation between the Kagome operators all 
and the pocket operators bU corresponds to a discrete 
Fourier transformation of a 2 x 2 lattice, and is given by 

ht3 — I Y" p^-aj f« - T. ia V (Q) 
17 Km — 2 / j a Km — / j M^Kmi \ v J 

I I 

where we have defined the symmetric and orthogonal 4x4 
matrix 

T . l=Tx . =T * l=T ji = le^a*. (10) 

Note that the matrix elements of T are ±1/2, as the 
scalar products Bj • a/ of Eq. (0 equal or ±7r. 

The tight-binding Hamiltonian is diagonal in the 
pocket indices j (cf. Appendix A Eq. (|A2)0 . 

#tb = £E e™ m 'bV m X m '„- (11) 

Ik. a mm' 

From this expression it is apparent, that the tight-binding 
Hamiltonian is invariant under any U (4) transformation 
of the of the form 

3' 

Eq. is just a special case of Eq. (|12|). This shows that 
i^tb is also diagonal in the Kagome indices. 

It is important to notice that the transformations in 
Eq. ([T2j) involves symmetries that are not present in a 
more general tight-binding model. For example a finite 
hopping integral tdd due to the a-hybridization between 
neighboring ^-orbitals would break this symmetry. We 
will discuss this aspect below in more detail and remain 
for the time being in this high-symmetry situation. 

In Na x Co02 the lower two bands are completely filled 
and will be quite inert. For this reason in the following 
sections we will only deal with the operators of the top 
band whose operators are denoted as 

«L = E°K m a^ and 6« =^0^ Kma , (13) 

m m 

resp., where 0^ m are matrix elements of the rotation 
matrix K of Eq. J3J). 

The top band gives rise to four identical Fermi surface 
pockets in the BZ, one in the T-point and three at the 
M-points. A translation in the reciprocal space by the 
vectors Bj maps the pocket around the Y point onto a 
pocket around the M-point. However, this fact does not 
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lead to nesting singularities in the susceptibility because 
a hole pocket is mapped onto a hole pocket by the vector 
Bj . The susceptibility of the top band is given by 



1 fk+q ~ fk 



Xq N ^ El 



N ^ Eh 



+ Q 



E 1 ' 

K + Q 



(14) 



where f k = f((3(El — /i)) and / is the Fermi-function. 
In the last expression of Eq. ([Hj). the sum over k is re- 
stricted to the reduced BZ. q also lies in the reduced BZ 
and is given by Q = q + Bj . The susceptibility Xq — Xq is 
periodic with respect to the reduced BZ and is just four 
times the susceptibility of a single Kagome lattice. As we 
have almost circular hole pockets with quadratic disper- 
sion around the T and the M points, the susceptibility is 
therefore approximately given by the susceptibility of the 
free electron gas in two dimensions within each reduced 
BZ, with circular plateaus of radius 2kf around the T 
and the three M points. 



III. COULOMB INTERACTION 

In this section we introduce the Coulomb interaction 
between the electrons. As we have spin and orbital de- 
grees of freedom, the on-site Coulomb interaction consists 
of intra-orbital repulsion U, inter-orbital repulsion U\ 
Hund's coupling Jh and a pair hopping term J'. These 
parameters are related by U ~ U' + 2 Jh and Jh = J 7 , 
where the first relation is exact for spherical symmetry. 
We can write the onsite Coulomb interaction as 



IK 



my^m' <j<j' 



9 C rrna C irn'a' Crm < j/Crrn ' (J 



(15) 



rriy^m' era' 



2 C rrna C \rna' C rrn'a'CYrn'(Ji 



where n rmcr = cj mcr c rmcr . We obtain an effective Hamil- 
tonian for the Coulomb interaction by rewriting the 
Hamiltonian in terms of the pocket operators of the top 
band defined in Eq. ([T3|). For small n = |k|<2 we can 
expand Eq. (|T3|) in powers of n 2 and obtain up to terms 
of the order h? 



12 



u K.mcr i 



(16) 



where m = 27rm/3. Expanding the energy of the top 
band around the point T', we obtain 

4 = t( 4 - J + ^ - ^ cos(60) + 0(* 8 )). (17) 

This shows that the pockets around the points V are al- 
most perfectly circular. The radius n^/a of these pockets 



9C = 


-3?7+2J / + 2Jh + 2?7 / 


9L> : 


= +3£/+6J / -2J H -2[/ / 


9£ c = 


+3?7-2J / -10Jh + 14[/ / 


9E S 


= -3[/+2J'-6J H + 2U' 


9F C = 


+3[/-2J / + 14J h -10[/ / 


9F S 


= -3f7+2J'+2J H -6*7' 



TABLE I: The coefficients of Eq. tfM 



depends on the Na doping x. Note that x corresponds 
to the density of carriers with x = 1 giving a completely 
filled top band. We have = 7r(l — x)/y/S. For the 
interaction in weak-coupling and at low temperatures, 
the states near the Fermi surface are important. For 
these states and for not too small Na doping x we can 
neglect the second term in the parenthesis of Eq. (|T6|) 
compared to 1. Note, that this condition on x is not 
very restrictive. Even for x = 0.35 the second term to- 
gether with all higher order terms is on the average one 
order of magnitude smaller than 1. Dropping the sec- 
ond term in Eq. ([TB|) spreads the a\ g symmetry of the 
states which is exact only for k = 0, to all relevant 
states in the top band. The interaction (p"5|) can now be 
rewritten in terms of the a\ g symmetric operators bH a . 
Processes involving states of the filled lower bands are 
dropped. The dropping of the second term in the paren- 
thesis of Eq. (p"6j) is a considerable simplification because 
it removes all K-dependence of the potential. 

At this point it is convenient to introduce density and 
spin density operators for the pocket operators of the top 
band. 



N ^ K 



+ QO" y KCT 



(18) 



The resulting effective interaction can be expressed with 
these operators in the following way 



N 



#eff = — ( B ijklSQS-Q 



^ijkl n Q n -Q 



(19) 



The symbols B c / S depend on the Coulomb integrals and 
are given by 



b: 



,c/s 
ijkl 



±C(25 ijk i - e ijkl ) ± D5uS jk 
+E c / s 5 ZJ 5 kl +F c / s 5 zk 5 jl , 



(20) 



where the S (e 2 ) symbol equals 1, if all the indices are 
equal (different) and otherwise. The coefficients C, D, 
E c / S and F c / S are listed in TABLE ffl Note, that for 
small pockets, the momenta k of the pocket operators 
&k in the four fermion terms of Eq (|TU|) can not add up 
to a half a reciprocal lattice vector b^. In order to con- 
serve momentum they must therefore add up to zero. 
Due to the position of the pockets in the BZ, Umklapp 
processes with low energy transfer are however possible 
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for arbitrary small pockets. In fact, the processes propor- 
tional to efj kl and 5u5jk(l — Sij) are Umklapp processes, 
as — Bj + &i — B/e is a non- vanishing reciprocal lat- 



U (4). The symbols A c J- kl defined in Eq. ((51)) are invariant 
under permutation of the indices, i.e. 



tice vector for e^/ ^ and for 5u5jk(l 



7^ 0, and 



aS/c 

Hjkl 



as/c 

^V{i)V{j)V{k)V(l) 



V eS 4 



(25) 



from Eq. © the momentum created by the operator bU From this follows that the symmetric group <S 4 is a sub- 



is K + Bj . 

Some details about the derivation of Eq. ((T§)) are pro- 
vided in Appendix El There are different ways of writing 
this interaction in terms of the operators in ((18)) . Our 
formulation treats charge- and spin degrees of freedom 
in the same way. It corresponds to the decomposition of 
a Hubbard interaction n^ni into \{\n 2 — S • S). 

In order to express the effective interaction Hamilto- 
nian of Eq. (JT9)l in terms of the Kagome operators a z Kcr , 
we define spin and charge density operators from the 
Kagome operators a l Kcr as in Eq. (|18|) . 



4 V- t* 



TV 

Kcr 

2 v-^ 



+ Q<7 U/ Kcr 



(21) 



Note, that the density operators, which are defined from 
the pocket operators b J Kcr are marked by a hut. The ef- 
fective Hamiltonian, i^ e ff? of Eq. (|19|) can be rewritten as 



H & = i E (^;«SgS-Q + ^. w ngn!_0 ■ (22) 



From Eq. (0) and JUH follows that 

^ijkl ~ Fim3~j n J-ko3~lpBj inop . (23) 

The symbols A c l s turn out to have a simpler structure, 
given by 



A c — — 



ijkl 



c 



(2J H - U f )S ik Sji 



(24) 



+ -^fiijkl — J'bilbjk — JuSijbkl — U'SikSji 



IV. SU(4) GENERATORS 

The tight-binding Hamiltonian described in section ITTI 
has a U (4) symmetry, reflecting the fact that it consists of 
4 independent and equivalent Kagome lattices. The cor- 
relations introduced by the on-site Coulomb repulsion in 
Eq. ((T5)) breaks this symmetry and leads to interaction 
between orbitals belonging to different Kagome lattices, 
as the three ^-orbitals on a given Co-site belong to three 
different Kagome lattices. The effective Hamiltonian in 
Eq. ((19)) is not invariant under general £7(4) transfor- 
mations, but is still invariant under a finite subgroup of 



group of G. Multiplying all operators a with the same 
Kagome index I by — 1 also leaves the Hamiltonian, i^ e ff, 
invariant, because the symbols A C J.^ are nonzero only if 
the four indices ijkl are pairwise equal. These two dif- 
ferent symmetry operations generate a group with 384 
elements. This group G is isomorphic to the symmetry 
group of the four-dimensional hypercube. In appendix O 
the structure of the group G is discussed and a character 
table is shown. 

To proceed, let Q r r = 0, . . . , 15 be a basis in the 16 
dimensional real vector space, V of Hermitian 4x4 matri- 
ces, fulfilling the usual ort honor mality and completeness 
relations 
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(26) 



r — 



This basis can be chosen such, that Q° is proportional 
to the unit matrix, (J 1-3 are diagonal, Q 4-9 are real and 
Q 10-15 are imaginary. It is convenient to define also the 
dual matrices 

-K-ij ~ 3~imJ- jnQmn- (27) 

In TABLE HH a choice of a basis Q r , which is particularly 
suitable for our purposes, is shown together with the dual 
basis K r . A representation p of the group G on V is given 
by p(g)Q r = NjQ r N g for ^GG, where N g is the natural 
four-dimensional representation of G (cf. appendix @ • 
The representation p is reducible and V is the direct sum 
of the four irreducible subspaces V°, V 1-3 , V 4-9 and 
V 10 - 15 spanned by matrices Q°, Q 1 " 3 , Q 4 ~ 9 and Q 10 " 15 , 
respectively. Therefore, the chosen Basis is appropriate 
for the symmetry group G. Defining charge and spin 
density operators 



n r Q = Q r ij n%=Kl j n% 

Sr _ /ar cij _ jy-r Aij 
Q — ^ij^Q — ^ij^Q 5 



(28) 



the interaction Hamiltonian can be written in a diagonal 
form as 



N 15 / 1 
H * = J E E ( A * S q S -q + 4 A > r Q n r 



r=0 Q 



(29) 



The coupling constants A^ s are equal for all Q r belong- 
ing to the same irreducible subspace in V. They are given 
in TABLE nm 



V. REDUCTION OF THE SYMMETRY 

The tight-binding Hamiltonian in Eq. (|11|) has a U (4) 
symmetry and even after introducing Coulomb interac- 
tion, the effective Hamiltonian ((25)) is invariant under the 
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TABLE II: The matrices Q 1 15 are a choice of an orthonormal complete basis of the 15 dimensional real vector space of traceless 
hermitian matrices, so called generators of SU(4), that is adequate to the symmetry of the CoCb-layer. The matrices K r are 
obtained from Q r by Eq. (l2?t. Note that 1 = -1 and i = -i. 2^2Q° = 2V2K is the 4 x 4 unit matrix. 



r 





1-3 


4-9 


10-15 




§(3*7+ W -6J H ) 


f(3[/-4[/ / + 2J H ) 


f(-2[/ / + 4J H -2J / ) 


§(-2[/ / + 4J H + 2J / ) 


A^ 


-|(3?7 + 6Jh) 


-|(3C7-2J h ) 


-%(2U' -2J') 


-f(2[/ / + 2J / ) 



TABLE III: The coefficients A^ 



symmetry group G. In a real C0O2 plane this symme- 
try is reduced even in the paramagnetic state. There are 
terms in the Hamiltonian of the real system that restrict 
the symmetry operations of G to the subgroup, which 
describes real crystallographic space-group symmetries. 



A trigonal distortion of the oxygen octahedra by ap- 
proaching the two O-layers to the Co-layer, is for exam- 
ple compatible with the point group symmetry Dsd of 
the C0O2 layer. However it lifts the degeneracy of the 
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i2g-orbitals, leading to a term 

H tx = D tl J2 J2 C kma C km'<r (30) 
ka m^m' 

— n t h\l hi , p iB Z -(a m -a m ,) 

Ik. (T m^m' 

in the Hamiltonian, where we used Eq. (|SJ to obtain the 
second line. For the top band we obtain with Eq. (fT5j) 
and (JB4|) 

H tr = v^A r 4^(^+0(^ 2 ))OL 

« y/2/3D tI Nni (31) 

where the matrix if 4 is given in TABLE ITTI and ft = |k|<2 
is small for the relevant states near the Fermi pockets, if 
the pockets are small enough. Similarly, a finite direct 
hopping integral tdd leads to the term 

H d d = tdd ^ 2 cos (fc m )4 m(7 Cfcmo- 
kma 

« V6t d(i 7Vn 4 (32) 

where we again dropped the terms involving the lower 
bands in the second line. In fact, any other additional 
hopping term or any quadratic perturbation compatible 
with the space group is proportional to the field n 4 in the 
limit of small pockets, if the perturbation is diagonal in 
the spin indices. As the trigonal distortion of the octahe- 
dra is nonzero and additional hopping terms are present 
in the Co02-layer, a term proportional to n 4 exists in the 
Hamiltonian acting like a symmetry breaking field. For 
simplicity, we will refer to a term proportional to n 4 in 
the Hamiltonian as the trigonal distortion, even though 
this term is rather an effective trigonal distortion that 
also includes the effects of additional hopping terms. 

From the matrix K A can be seen, that the presence of 
a finite field, n 4 , in the Hamiltonian leads to a distinc- 
tion between the Y and the M points in the BZ and the 
four hole pockets are no longer equivalent. In real space, 
the four Kagome lattices are still equivalent, as they 
transform under space group symmetries among them- 
selves. In fact, the matrix Q 4 is still invariant under 
permutations of rows and columns, i.e. N^Q 4 N g = Q 4 
for all g G £4, but Q 4 is not invariant under changing 
the sign of all operators with the same Kagome index. 
These sign changes, however, are not space-group sym- 
metries, but gauge symmetries, originating from the fact 
that the charge on the Kagome lattices is conserved by 
Htb and also by the Coulomb interaction except for the 
pair-hopping term proportional to J' in Eq. ([T5J). This 
term however can only change the number of electrons by 
two, leading to these gauge symmetries, that are broken, 
as soon as single electron hopping processes between the 
Kagome lattices are introduced. 



To classify the states according to the real symmetry 
group of the Co02-layer without gauge symmetries, it is 
therefore sufficient, to consider the presence of a small 
field n 4 , that restricts the symmetry group G to a sub- 
group, consisting of space group symmetries of the C0O2- 
layer. This subgroup of G is isomorphic to S4 c± Td — O. 
Intuitively it is understandable that the symmetry of the 
four dimensional cube reduces to the symmetry of a three 
dimensional cube, if one of the four hole pockets is not 
equivalent to the other three. 

Form TABLE ITTI can be seen, that the matrices Q°, 
Q 1 - 3 , Q 4 , Q 5 - 6 , Q 7 - 9 , Q 10 - 12 and Q 13 ~ 15 transform ir- 
reducibly under £4 with the representations TJ, Tg, 
T3, r|, T4 and r§, respectively, where the upper-script 
letter distinguishes between different subspaces trans- 
forming with the same representation. 

The appearance of three dimensional irreducible rep- 
resentations in the classification of the order parameters 
can be understood as follows. The point group P of a sin- 
gle Co02-layer is Dsd^ and the degree of its irreducible 
representations is < 2. The point group is the factor 
group S/T where S is the space group of the C0O2- 
layer and T is the subgroup of all pure translations. For 
our system it is convenient to consider the factor group 
P' = S/2T, where 2T is the subgroup of T that is gener- 
ated by translations of 2a^. P' is isomorphic to the cubic 
group Oh and has irreducible representations of degree 3. 
The operators and Sq transforms irreducibly under 
the translations in 2T for every r. The symmetry opera- 
tions of P' however mix operators (or Sq) with differ- 
ent r, and the irreducible representations as given above 
or shown in TABLE [H] are obtained. Strictly speaking, 
the basis of SU(4) generators shown in TABLE ITTI is the 
correct eigenbasis only for an infinitesimal small trigonal 
distortion, for a finite distortion, the representations Y\ 
and Y\ as well as and Y\ can hybridize as they trans- 
form with the same irreducible representation. Note, that 
r§ transforms differently under time reversal. The situ- 
ation here is similar to atomic physics, where a crossover 
from the Zeeman effect to the Paschen-Back effect with 
increasing magnetic field occurs, because states with the 
same J z can hybridize. 



VI. ORDERING PATTERNS 

In this section the different types of symmetry break- 
ing phase transitions are discussed in a mean-field pic- 
ture. The symmetry breaking is due to existence of a 
finite order parameter, that is in our case given by the 
expectation value (n^) of (Sq). Note, that a finite expec- 
tation value (tIq) or (n 4 ) does not break any symmetry 
of the Co02-layer. 

In our tight-binding model as it was discussed in sec- 
tion |H] the susceptibility, x° is given by 4 identical 
plateaux around the T and the M points. In the presence 
of a trigonal distortion, the susceptibility still keeps a 
plateaux like structure but the diameter of the plateaux 
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decreases, such that the susceptibility appears sharply 
enhanced around the M and the T points. Therefore we 
restrict the discussion to the case where Q equals zero 
and write n r and S r instead of n r and S£ from now on. 
Note, that in our formalism the states with q = de- 
scribe periodic states with the enlarged unitcell of the 
Kagome lattice. But the internal degrees of freedom 
within this enlarged unitcell still allows for rather com- 
plicated charge- and spin-patterns. States with a small 
but finite Q describe modulations of these local states 
on long wavelengths. It is therefore important to under- 
stand first the local states the are described by Q = 
instabilities. Furthermore, only q = states couple to 
the periodic potential produced by a Na superstructure 
at x = 0.5. 

The Q = instabilities lead to a chemical potential dif- 
ference for states belonging to different hole pockets. In 
general, the BZ is folded and states of different hole pock- 
ets combine to new quasi-particles. In this case, trans- 
lational and/or rotational symmetry is broken. Com- 
plex ordering patterns can be realized without opening 
of gaps, i.e. the system stays metallic. 

We consider first the orderings given by a finite ex- 
pectation value of the charge density operators n r . This 
expectation value is given by 

K) = |E A ^k^L>, (33) 

Kcr 

where A[ are the eigenvalues of the matrix Q r 
(UkiQijUfj = A[fc) and v l Ka = Uf n a™ a are the creation 
operators of the quasi-particles. If only one (n r ) ^ 0, 
the effective interaction Hamiltonian in the mean-field 
approximation reduces to 

^E A F^^- (34) 

Kit 

If the coupling constant K Q r is negative, the interac- 
tion energy of the system can be lowered by introduc- 
ing an imbalance between the occupation numbers rt\ = 
^Kcr( v Kcr^Lcr)- The operators v\^ s create Bloch states 
with momentum k in the reduced BZ. The amplitudes of 
the three ti g orbitals on a given Co site with these Bloch 
states can be obtained from Eq. (0) and {7J) and the rela- 
tion <2k ~ 1/ a/3 CL^m which follows from Eq. ® and 

m- 

For the matrices Q° 4 , these Bloch states are given by 
a single ti g orbital on each Co site. For the non-diagonal 
matrices Q 4 ~ 9 these Bloch states are on each Co site 
proportional to a linear combination of t^g orbitals of 
the form 

A=(s x d x + Sydy + s z d z ) with s x , s y , s z e {±1}. (35) 
V3 

This linear combination is the atomic <i-orbital cpo = J^o 
parallel to the body-diagonal [s x , s yi s z ] of the cubic unit- 
cell around a Co atom. 



The eigenvectors of the of the matrices Q 10-15 are com- 
plex. A complex linear combination of t<i g orbitals has in 
general a non-vanishing expectation value of the orbital 
angular momentum operator L. In TABLE IIVI the an- 
gular momentum expectation values, which are relevant 
for our discussion are shown. 



(id x + d y + d z )/y/3 


(L) 


= ^(0,-1,1)2/3 


(cyclic) 


(id x + d y - d z )/V3 


(L> 


= 71(0,1,1)2/3 


(cyclic) 


(d x + uj 2 d y + cud z )/V3 


<L> 


= i, i)/V3 




(ou 2 d y +(jud z )/V2 


(L> 


= fc(l,0,0)>/3/2 


(cyclic) 



TABLE IV: The expectation values the angular momentum 
operator L for several complex linear combinations of t2 g or- 
bitals. ou = e 2irt/s 

The quasi-particles v l K(T are expressed in terms of 
pocket operators by v l Ka = Ui m b™ a , where the unitary 
matrix U[ m = Ul n T n m diagonalizes K r . From this fol- 
lows that if K r is already diagonal, no folding of the BZ 
occurs and translational symmetry is not broken. Oth- 
erwise, the BZ is folded and states of different pockets 
recombine to form the new quasi-particles. 

Now we consider finite expectation values of the spin- 
density operators S r . Due to the absence of spin-orbit 
coupling, our model has an SU(2) rotational symmetry 
in spin space. Therefore the discussion can be restricted 
to the order parameters (SI) = (e z • S r ), given by 

^) = |E A ^(^L), (36) 

Kcr 

where a takes the values 1 and —1 corresponding to spin 
up and down. If only one ^ 0, the effective interac- 
tion Hamiltonian reduces to 

— 2 — Z^ x i cr v L* v k* • (37) 

K(J 

The mean-field Hamiltonian (pTfj) is given by the same 
quasi-particles and the same eigenvalues A[ as the Hamil- 
tonian in (|34[) . The only difference is that the sign of 
the splitting of the quasi-particle bands depends on the 
spin. In the following, all ordering transitions with order 
parameters (n r ) and (SI) for r = 0, . . . , 15 are shortly 
discussed. 



r=0 

charge: (n°) is the total charge of the system, which is 
fixed and non-zero, even in the paramagnetic phase. 

spin: A finite (S®) describes a Stoner ferromagnetic in- 
stability. The coupling constant Aq given in TABLE IlIII 
is the most negative coupling constant. In the unper- 
turbed system without trigonal distortion, the critical 
temperature of all continuous transitions discussed here, 
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only depends on the density of states and on the cou- 
pling constant in the mean-field picture. In this case, 
ferromagnetism is the leading instability for the unper- 
turbed system. In the real Co02-plane, this must not 
necessarily occur, but strong ferromagnetic fluctuations 
will be present in any case. 



r=l 3 

charge: A finite expectation value (rf) for r = 1,2,3 
corresponds to a difference in the charge density on the 
four Kagome lattices, because the matrices Q 1-3 of TA- 
BLE [D] are diagonal and the quasi-particles v^ a are just 
the Kagome states aUa- From the view point of Fermi 
surface pockets given by if 1-3 which are non-diagonal, 
this order yields a folding of the BZ, because the quasi- 
particles vlla- are linear combinations of states belonging 
to different hole pockets. This means that the trans- 
lational symmetry is broken. In the matrix Q 1-3 we 
find two positive and two negative diagonal elements. 
Consequently, a finite expectation value (n 1-3 ) leads 
to a charge enhancement on two Kagome lattices and 
to a charge reduction on the other two. As specifying 
two Kagome lattices specifies a direction on the trian- 
gular lattice, rotational symmetry is broken and crystal 
symmetry is reduced from hexagonal to orthorhombic. 
The phases described by the matrices Q 1-3 have the 
same coupling constant A^ because they transform ir- 
reducibly into each other under crystal symmetries with 
the representation Tg. In order to examine which linear 
combinations of the three order parameters (n 1 ), (n 2 ) 
and (n 3 ) could be stable below the critical temperature, 
we consider the Landau expansion of the free energy 
AF = F - F 



AF 



72 



4 



(vW 2 



4 



% 2 ) 2 

(38) 



with rji = (n 1 ), 772 = ( n<2 ), V3 — (tt 3 ). For 71 > 
max{0, — 72}, the free energy is globally stable. For 
72 < 0, Eq. (|55J) has a minimum of the form 771 = 772 = 773, 
if /3 2 — 4a (371 +72) > 0. This phase is described by the 
symmetric combination Q 1 = (Q 1 + Q 2 + Q 3 )/V3 which 
does not break the rotational symmetry. In FIG. 03 the 
folding of BZ and the splitting of the bands (the dotted 
line is triply degenerate) and the orbital pattern of the 
quasi-particles vll a = a^ a are shown. Note that Q 1 has 
one positive and three negative diagonal elements. The 
charge is enhanced or reduced on a single Kagome lattice 
depending on the sign of the coefficient f3 in Eq. The 
third-order term in the free energy expansion is allowed 
by symmetry, because there is no inversion- like symmetry 
that would switch the of (771,772,773) — > (—771, —772, —773). 
Therefore the transition can be first order. On the other 
hand, for 72 > 0, there is a competition between the 
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FIG. 3: Charge ordering instability with finite expectation 
value (ho) leading to a charge enhancement or reduction on 
one Kagome lattice. The folding of the BZ and the splitting 
of the pockets is shown (The double dotted line in the BZ 
indicates a triply degenerate pocket). On the right a quasi- 
particle state that is in this case just a Kagome lattice state 
is drawn. 



terms proportional to 72 and (3 in Eq. (|58J). The mini- 
mum has not a simple form. For \/3\ <C 72, however, the 
transition yields states approximately described by the 
matrix Q 1 , Q 2 or Q 3 . In any case this phase does break 
the rotational symmetry. 

spin: The spin density mean-fields (S l z ) i = 1,2,3 
transform under space group symmetries like and time 
reversal symmetry gives (S l z ) to —(SI). Due to the latter 
the third order term in Eq. ((38)) is forbidden, so that the 
transition is continuous. For 72 < 0, Eq. (|58J) has again 
a minimum of the form 771 = 772 = 773, whereas for 72 > 
the minimum is realized for 771 7^ and 772 = 773 = (and 
permutations) , if a < 0. The folding of the BZ, the 
quasi-particles and the breaking of space-group symme- 
tries is the same as for the charge density operators n\ 
However, the splitting of the bands depends now on the 
spin and time reversal symmetry is broken. 

These states are spin density waves, spatial modula- 
tions of the spin density with a vanishing total magnetiza- 
tion. The two different types of spin density modulations 
for 72 > or 72 < are shown in Fig. 21 For 72 > ro- 
tational and translational symmetry is broken yielding a 
collinear spin orientation along one spatial direction and 
alternation perpendicular. In contrast 72 > yields a 
rotationally symmetric spin density wave with a doubled 
unit cell. This special type of spin density wave gives 
a subset of lattice points, forming a triangular lattice, 
of large spin density and another subset with opposite 
spin density of a third in size, forming a Kagome lattice. 
Both states are metallic, because no gaps are opened at 
the FS. This spin density wave is not a result of Fermi 
surface nesting, but due to the complex orbital struc- 
ture. The coupling constant for this transition, A\ is the 
second strongest coupling in the model Hamiltonian af- 
ter the ferromagnetic coupling constant, Aq, as it is best 
seen in Fig. 
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FIG. 4: The spin density wave patterns corresponding to a 
finite expectation value (Si) is shown on the left. This pattern 
is stabilized if 72 > in Eq. (|38|>. The pattern on the right 
corresponds to a finite order parameter (Si + + S^), which 
is stabilized for 72 < 0. 



r=4 



charge: As discussed in section a finite expectation 
value of n 4 does not break any space group symmetry. 
The matrix K A is diagonal with one positive and three 
negative elements. This leads to a change of the band 
energy of the band at the Y point relative to those at 
the M points (Fig. Et-))- This results in an orbital or- 
der, a pattern as shown in Fig. Ol), because the number 
of holes associated with the hole pocket around the Y- 
point is different from that of the other pockets. The net 
charge onsite vanishes, but the charge distribution has 
the quadrupolar form, which results from 



p(r) oc j [3 1 ^ + ipzx + ^ xy \ 2 ~ 

~\^yz +^zx ~^xy\ 2 ~ |VV' 
= ^y Z ^zx + Iptxlpxy + ^ty^yx 



\lpyz ~ Ipzx - 

:-^ZX~ Ipx 



• ^X 



c.c. 



(39) 

The corresponding tensor operator belongs to the repre- 
sentation Yi of the subgroup D3 of the cubic group with 
the three-fold rotation axis parallel to [111], i.e. along 
the c-axis perpendicular to layer. This quadrupolar field 
would be driven by the symmetry reduction discussed 
above, through trigonal distortion and direct (id-hopping 
among the ^-orbitals. 

spin: While the corresponding order parameter (Sf) 
breaks time reversal symmetry, space group symmetry 
is conserved. This order is spatially uniform analogous 
to a ferromagnet without, however, having a net mag- 
netic moment. Because the magnetic moment associated 
with the Fermi surface pocket at the T-point is oppo- 
site and three times larger than the moment at the three 
M-pockets. While the net dipole moment vanishes on ev- 
ery site, this configuration has a finite quadrupolar spin 
density corresponding to the onsite spin density distri- 
bution of the same form as the charge distribution in 
Ea. ((59)1 . which also belongs to Yi representation of D 3 . 
It is also important to note that no third order terms are 
allowed due to broken time reversal symmetry, such that 
the transition to this order would be continuous. 




FIG. 5: Ordering instabilities, described by real off-diagonal 
SU(4) generators Q 4-9 . a) Q 4 breaks neither translational 
nor rotational symmetry, b) shows the ordering corresponding 
to Q 6 . The ordering shown in c) is described in real space and 
in reciprocal space by the same matrix Q 7 = K 7 and breaks 
translational symmetry. The corresponding BZ is shown in 
FIG. El 



r=5,6 

charge: The order parameters (n 5 ) = 775 and (n 6 ) = t]q 
transform according to the irreducible representation T3 
of the cubic point group. The Landau expansion of the 
free energy is given by 



AF 



a 



■nl) 



3% (3% 



•% 2 ) + 7(^ + %T> (40) 



whose global stability requires 7 > 0. The third or- 
der term, allowed here, induces a first order transition 
and simultaneously introduce an anisotropy which is not 
present in the second- and fourth-order terms. We can 
write (775 , 776) = ^(cos tp, sin (p) and obtain 



AF 



-77 + -77 sin3<p+ -77 



(41) 



Depending on the sign of j3 the stable angles will be ip = 
sign(/?)7r/2 + 27rn/3. This yields three degenerate states 
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of uniform orbital order whose charge distribution has 
the quadrupolar form: 

p(r) oc e ltp {{^* zx ^ xy + iplyipzx) + u{$* yz $ zx + ipl x ipyz) 

(42) 

with a tensor operator belonging to T3 of D3. Each state 
is connected with the choice of one M-pocket which has 
a different filling compared to the other two (Fig. Eb). 
The main axis of each state points locally along one of the 
three cubic body-diagonals, [1, 1, 1], [1, 1, 1], [1, 1, 1], and 
the sign of the local orbital wave function is staggered 
along the corresponding direction on the triangular lat- 
tice, [2, 1, 1], [1, 2, 1], [1, 1, 2). In this way the rotational 
symmetry is broken but the translational symmetry is 
conserved. The matrices Q 5 and Q 6 commute with Q 4 
such that the external symmetry reduction has only a 
small effect on this type of order. 

spin: The spin densities (5f ) and (S®) also belong to 
the two-dimensional representation Ts of the cubic point 
group. Here time reversal symmetry ensures that the 
Landau expansion only allows even orders of the order pa- 
rameter (775,776) = 7?(cos(/?, sin<^) (i.e. (3 = in Ea. ((l0j) . 
The continuous degeneracy in cp is only lifted by the sixth 
order term, given by 



2\ 2 w J- t> 1 w z 'I • 'I o 



2 ^2 . 



6 



6 



(43) 

Stability requires Si > max{0, —£2}. The anisotropy is 
lifted by the ^-term which give rise to two possible sets 
of three-fold degenerate states. Depending on the sign 
of 82 we have a minimum of the free energy for ip — 
(1 — signed )7r/4 + irn. The corresponding spin densities 
have no net dipole on every site, but again a quadrupolar 
form of the same symmetry as for the charge, given by 



r=7-9 
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FIG. 6: Transitions to time reversal symmetry breaking 
states, where the expectation value of the orbital angular mo- 
mentum (L) on the Co-sites is finite, a) States where the 
angular momentum does not lie in the plane. K 10 — Q 10 b) 



States with angular momentum in the plane. K 



The folding of the BZ and the hybridization of the bands for 
a) are shown. Dotted lines indicate double degenerate bands. 



charge: The order parameters (n l ) for i = 7, 8, 9 trans- 
form irreducibly under space group symmetries with the 
representation r|. The expansion ([58)1 of the free en- 
ergy holds also for these order parameters. The third 
order term makes the transition first order and favors 
the symmetric rotationally invariant combination of the 
order parameters, described by the matrix Q 7 = (Q 7 + 
Q 8 + Q 9 )/V3 = K 7 shown in FIG. Et). The folding of 
the BZ and the splitting of the bands is the same as in 
Fig. 03 The orbital pattern of the non-degenerate quasi- 
particle band is also shown in Fig. Et). It consists of 
atomic <p° orbitals pointing along all four cubic space di- 
agonals. Translational but not rotational symmetry is 
broken. 

spin: The discussion for the spin density operators is 
analogous to the discussion in the section r = 1 — 3. 



r=10-12 

charge: The order parameters (n l ) for i = 10,11,12 
transform irreducibly under space group symmetries with 
the representation T4. For the T4 representation of X^, 
there is no third order invariant. All other terms in 
Eq. (|38|) are however also invariants for T^. The absence 
of the third order term leads to continuous transition. 
The stabilized state for a < depends on the sign of 72 
in Eq. (158)) . 

For 72 > a nontrivial minimum with (n 11 ) = (n 12 ) = 
exists, which is described by the hermitian, imagi- 
nary matrix Q 10 . If A is an eigenvalue of Q 10 , then 
—A is also an eigenvalue of Q 10 and the corresponding 
quasi-particles are connected by time reversal symmetry. 
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Therefore the non- vanishing eigenvalues of Q 10 belong to 
quasi-particle states, which are not invariant under time 
reversal symmetry. They are given by complex linear 
combinations of ti g orbit als. For complex linear combi- 
nations of t2 g orbitals, the expectation value of the orbital 
angular momentum operator (L) does not vanish in gen- 
eral, as can be seen from TABLE HVl In FIG.|HK)? the 
pattern of the angular momentum expectation values (L) 
for a quasi-particle of Q 10 is shown. It is invariant under 
translations along a x and staggered under translations 
along a2 and a3. The expectation values are parallel to 
[Oil]. The folding of the BZ and the splitting of the 
bands is shown in Fig. Eh). Rotational, translational and 
time reversal symmetry is broken and the state has the 
magnetic point-group 2/ra. 

For 72 < the symmetric combination Q 10 = (Q 10 + 
Q 11 + Q 12 )/V3 = K 10 is stabilized. The angular mo- 
mentum pattern for a quasi-particle with non-vanishing 
eigenvalue is shown in FIG.[5Ji). Depending on the site, 
the expectation value points along [100], [010], [001] or 
[111] and the magnitudes are such, that the pattern is ro- 
tationally invariant and the expectation value of the total 
angular momentum perpendicular to the plane vanishes. 
The folding of the BZ and the splitting of the pockets 
is shown in FIG. Eh). This state has the magnetic point 
group 3m. Note, that these states can also be considered 
as a kind of staggered flux states. The matrices Q 10-12 
commute with Q A and therefore the transitions are only 
little affected by a trigonal distortion. 

spin: The spin density order parameters (S l z ) for i = 
10, 11, 12 also transform under space group symmetries 
like T 4 and except for the spin dependent quasi-particle 
energy, the discussion is the same as for the charge den- 
sity operators. Note, however, that these spin density 
operators do not change sign under time reversal sym- 
metry, because both the orbital angular momentum and 
the spin is reversed. This, however does not lead to a 
third order term in the Landau expansion, as there is no 
third order invariant for the T4 representation anyway. 



r=13-15 

charge: The order parameters (ri 1 ) for i = 13, 14, 15 
transform irreducibly under space group symmetries with 
the representation r§. The matrices Q 13-15 are also 
imaginary and time reversal symmetry changes the sign 
of the order parameters. The Landau expansion of the 
free energy is given as above by Eq. with (3 = 0. 

For 72 > and a < a minimum of the free energy is 
given by the order parameter (n 13 ). The angular momen- 
tum pattern of the quasi-particles is shown in FIG.Et>). 
The expectation values lie in the Co02-plane and are 
parallel to the ai direction. Their sign is staggered along 
the a2 and a3 direction. The quasi-particles consist of 
states belonging to the T and the M pocket. The folding 
of the BZ is given in FIG. [St), but with the single dotted 
line in the center being a doubly degenerate M pocket. 



Rotational, translational and time reversal symmetry is 
broken. 

For 72 < the symmetric combination Q 13 = (Q 13 + 
Q 14 + Q 15 )/V3 = -K 13 is stabilized. The pattern 
of the quasi-particles corresponding to Q 13 is shown in 
FIG. Eh). It consists of non-magnetic sites with a (po or- 
bital perpendicular to the plane and of sites with angular 
momentum expectation values along a^. Rotational sym- 
metry is not broken in this case. The folding of the BZ 
and the splitting of the bands is shown in FIG.Efc). 

All angular momentum expectation values for these 
two states lie in the Co02-plane. Therefore, it is not 
possible to interpret these states as staggered flux states. 

spin: The spin density order parameters (S l z ) for 
i = 13, 14, 15 are invariant under time reversal symmetry. 
Therefore, the third order term in Eq. (|38|) is allowed and 
the transition is a first order transition. 



VII. POSSIBLE INSTABILITIES 
A. coupling constants 

As can be seen from TABLE IIIII the coupling con- 
stants for the SDW transitions A s are rather negative 
whereas the charge coupling constants A c tend to be pos- 
itive. This is not surprising as only local repulsive inter- 
action is considered here, that tends to spread out the 
charge as much as possible. 

The coupling constants A^ s with r = 0, . . . , 3 depend 
on the intra-orbital Coulomb repulsion U . As U is the 
largest Coulomb integral, the absolute value of these cou- 
pling constants is biggest. The remaining coupling con- 
stants A^ s do not depend on U. For J' = they are also 
independent of r. For finite J' the degeneracy between 
the real (4-9) and imaginary (10-15) SU(4) generators is 
lifted. 

In order to compare the different coupling constants 
better, the relations U = U' + 2 Jh and Jh = J', that 
hold in spherically symmetric system, can be assumed 
to hold approximately. The ratio a = U' /U is positive 
and usually larger than 1/2 and smaller than 1. These 
assumptions allow to order the dimensionless coupling 
constants A c / S = 9A c / s /(2L r ) according to their strength. 
In FIG. [3 the dimensionless coupling constants A£ are 
shown as functions of a. The most negative coupling con- 
stant is the ferromagnetic one with Aq = — 6 + 3a. For 
a close to 1, the coupling constant for spin density order 
Af = — (2 + a) is comparable. Smaller but still clearly 
negative are also the coupling constants for the spin den- 
sity angular momentum states A\ = — (1 + a). The 
coupling constants A| = A| = 1 — 3a are also negative. 
Finally, the coupling constant for time reversal symmetry 
breaking angular momentum states A^q = 3 — 5a and for 
the charge density order A^ = 4 — 5a are rather positive, 
but can in principle also be negative if a is close enough 
to one. In fact it is quite remarkable that for a > 0.8 
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(0,s) 

a=0.5 a=1.0 

FIG. 7: The dimensionless coupling constants A^ s = 
9A c/ 7(2[/) as functions of a = U'/U. The relations U = 
U' +2 Jh and J' = Jh are assumed to hold. The solid (dashed) 
lines denote the charge (spin) coupling constants. 

all coupling constants constants (except A§)are negative. 
For a = 1 additional degeneracies among the coupling 
constants appear, as can be seen in FIG. This indi- 
cates the existence of a higher symmetry at this point. 
In fact, the local Coulomb interaction Hf? of Eq. [El de- 
pends only on the total charge n r = ^2 mcr n rmcr on the 
site r and is given by Un r (n r — l)/2 for a = 1. 



are finite, would be strongly affected by the trigonal dis- 
tortion, as the pocket states that hybridize in such a tran- 
sition are no longer degenerate. 

Finally, as mentioned above, the order parameters de- 
scribed by the matrices if 1-3 and K 7-9 transform with 
the same representation and are mixed by the trigonal 
distortion. For strong distortions the mixing tends to 
odd-even combinations and only the odd combinations, 
K 1 - K 7 , K 2 - if 8 , K 3 - K 9 commute with the symme- 
try breaking field, if 4 , and connects the still degenerate 
states of the M pockets. 

If the trigonal distortion is so strong, that the pockets 
states at the M points lie below the FS, only a sponta- 
neous ferromagnetic instability can still occur according 
to the Stoner criterion. First order transitions, however, 
are still possible. 

The ferromagnetism is the leading instability in the 
symmetric model and is least affected by the trigonal 
distortion. Therefore, in real Na^Co02 systems where a 
rather strong trigonal distortion is unavoidably present, 
ferromagnetism would be most robust and is in fact the 
only type of all the described, exotic symmetry breaking 
states, that would have a chance to occur spontaneously. 

However, even if the coupling constants of the more ex- 
otic states are not negative enough, to produce a sponta- 
neous instability, their corresponding susceptibilities can 
be large enough to give rise to an important response of 
the electrons in Co02-plane to external perturbations. 
In the next section, we describe how te Na-ions can be 
viewed as an external field for the charge degrees of free- 
dom. 



B. effect of the trigonal distortion 

In the mean-field description, an instability occurs if 
the Stoner-type criterion is satisfied. At zero tempera- 
ture in the system with full symmetry, this criterion reads 
in our notation as 

-A c/s 

—f-DiE*) = 1, (44) 

where D(E-p) is the density of states per spin and per 
hole pocket. For rather small pockets D(E-p) is given by 
y/3/(27rt) « 0.28/t in our tight binding model, increases 
however with decreasing Ep (cf. Fig. EJ). From Eq. 
we can estimate that the critical U must be larger than 
lOt for having a ferromagnetic instability. With the in- 
troduction of the trigonal distortion, as it was discussed 
in section El the Stoner criteria of Eq. (@2J) are modified. 

For the order parameters described by the matrices 
K°, K A ,K b ~ Q and if 10-12 , that commute with the trig- 
onal distortion if 4 , the change of the Stoner criterion is 
only due to the changing of the density of states at the 
M and the T pockets by the trigonal distortion, and the 
Stoner criterion is only slightly modified as long as all 
four pockets exist. 

On the other hand, the instabilities towards states, 
where the order parameters with the matrices K 13-15 



VIII. NA-SUPERSTRUCTURES 

In Na^Co02 the Na-ions separate the C0O2 planes. 
There are two different Na-positions which are both in 
prismatic coordination with the nearest O-ions. The Na2 
position is also in prismatic coordination with the near- 
est Co-ions, while the Nal position lies along the c-axes 
between two Co- ions below and above. This leads to 
significant Na-Co repulsion, suggesting that the Nal po- 
sition is higher in energy. In fact, the Na2 position is the 
preferred site for Nao.75Co02, where the ratio of occupied 
Nal-sites to occupied Na2-sites is about \\2w2k Deinter- 
calation of Na does however not lead to a further deple- 
tion of the Nal-sites. On the contrary, the occupancy 
ratio goes to 1 for x going to 0.5. Further there is a 
clear experimental evidence, that at x — 0.5 the Na-ions 
form a commensurable orthorhombic superstructure al- 
ready at room temperature^ For several other values of 
x also superstructure formation has been reported, but 
x = 0.5 shows the strongest signals and has the sim- 
plest superstructure Aifi In addition for x = 0.5 samples 
a sharp increase of the resistivity at 50K respectively at 
30K was reported iAL2i 

This experimental situation is rather surprising, 
naively one expects commensurability effects to be 
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FIG. 8: Two different Na superstructures in Na^CoC^- The 
left one does not break rotational symmetry and would drive 
a charge ordering as shown in FIG. The right one is in 
fact realized in Na^CoCb, it is obtained from the right one 
by shifting the Na-chains along the arrows. This shift is due 
to the Coulomb repulsion of the Na-ions. 




FIG. 9: The charge ordering pattern corresponding to the 
matrix K 1 —K 7 , consists of alternating rows of d x and d y — d z 
orbitals. On the left hand side the original BZ, the orthorhom- 
bic BZ due to the charge ordering and the experimentally 
observed reduced orthorhombic BZ (dark) are shown. 



strongest at x = 1/3 or at £ = 2/3 on a triangular lattice 
but not at x = 1/2. Therefore, it was concluded that 
structural and electronic degrees of freedom are coupled 
in a subtle manner in Na^CoC^™ 

In this section we show how the different ordering pat- 
terns can couple to the observed Na superstructure at 
x = 1/2. Before going into the details, we note that 
due to our starting point of inter-penetrating Kagome 
lattices, commensurability effects will be strongest for 
samples where the Na-ions can form a simple periodic 
superstructures that double or quadruple the area of the 
unitcell, since specifying a single Kagome lattice also 
quadruples the unitcell. For x = 1/2 such simple su- 
perstructures exist as shown in FIG. [HI A sodium super- 
structure couples to the charge but not to the spin de- 
grees of freedom in the C0O2 layer. In our model, there 
are 15 collective charge degrees of freedom. From FIG.[7| 
can be seen that A£ is most negative for r = 4, . . . , 9. 
Hence, these modes are the "softest" charge modes gen- 
erating the strongest response to a Na-pattern. As shown 
in FIG. [SJ the charge order corresponding to r = 4, 5,6 
does not enlarge the unitcell and does therefore not opti- 
mally couple to the Na-patterns that can be formed with 



x = 0.5. However the orbital pattern shown in FIG.JSt) 
has lobes of electron density pointing towards selected 
Nal and Na2 positions. For x — 1/2 it is possible to 
occupy all these and only these positions. This leads to 
the left Na-superstructure of FIG. [HI In other words, this 
Na-superstructure couples in a optimal way to this rota- 
tionally symmetric charge pattern. Further, the Landau 
expansion shows that the rotationally symmetric combi- 
nation is favored by the third-order term. Therefore, it 
is clear that the electronic degrees of freedom would fa- 
vor this Na-superstructure. This pattern however does 
not maximize the Na-ions distances. It is apparent that 
the average distances between the sodium ions can be in- 
creased, if every second of the one-dimensional sodium 
chains is shifted by one lattice constant as shown in 
FIG. [HI In this way an orthorhombic Na-superstructure 
is obtained, which is the one observed in experiments. 
This orthorhombic pattern does not drive the rotation- 
ally symmetric charge pattern shown in FIG.[5t), which 
is described by the matrix K 7 = K 7 + K 8 + if 9 , it might 
however drive the orthorhombic charge pattern described 
by the matrix K 7 or rather the orthorhombic charge pat- 
tern described by K 1 — K 7 , as in the presence of trigonal 
distortion the K 1 and K 7 mix and the odd combina- 
tion will have the most negative coupling constant. This 
charge pattern is shown in FIG. 03 It consists of lines 
of c^-orbitals alternating with lines of the linear combi- 
nation d y — d z -orbitals. Note, that this charge pattern 
corresponds to the mixed K 1 — K 7 -matrix, the charge 
is not uniformly distributed on the Co-atoms. In this 
charge pattern, the Nal-sites above the (d y — d z ) Co- 
sites will be lower in energy than the Nal-sites above the 
d x Co-sites and similarly the Na2-positions are separated 
into nonequivalent rows. 

In reciprocal space, such a charge ordering leads to 
a folding of the BZ such that the two M-pockets hy- 
bridize. The ordering of the Na-ions along the chains 
leads to a further folding of the BZ and to a hybridiza- 
tion of the bands, as it is shown in FIG. [51 The schematic 
FS in FIG. El is drawn to illustrate the hybridization 
occurring due to the translational symmetry breaking. 
Li et al. performed density-functional calculations in or- 
der to determine the band-structure of Nao.sCoC^ in the 
presence of the orthorhombic superstructure from first- 
principles. 49 Quite generally one can assume that this 
superstructure, which specifies a direction on the trian- 
gular lattice, can lead to quasi-onedimensional bands in 
the reduced BZ. For such one-dimensional bands, nesting 
features are likely to occur and would lead to a SDW-like 
instability, as it was observed at 53 K by Huang et al.£2£ 
Such a transition could open a gap at least on parts of 
the FS and in this way lead to the drastic increase of 
the resistivity observed at 53 Ki At higher temperature, 
the resistivity is comparable in magnitude to the metallic 
samples and increases only slightly with lowering temper- 
ature. This weakly insulating behavior could be another 
effect of Na-ion ordering. Since the rotational symme- 
try is broken, domains can be formed. The existence of 
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domain walls would be an obstacle for transport where 
thermally activated tunneling processes play a role. It 
would be interesting to test this idea by removing the do- 
mains and see whether metallic temperature dependence 
of the resistivity would result. A bias on the domains can 
be given by in-plane uniaxial distortion. 

To finish this section, we will discuss a further mecha- 
nism, that could lead to a non-magnetic low-temperature 
instability in Nao.sCoC^. In section IVII we saw that 
the third order term in the Landau expansion, Eq. 
favors always a rotationally symmetric charge ordering 
where all three order parameters 771, 772 and 773 have the 
same magnitude. But as argued above, the Na-ion repul- 
sion leads nevertheless to an orthorhombic charge order- 
ing, where only one order parameter 771 is finite. From 
Eq. we obtain a Landau expansion for the remaining 
two order parameters 772 and 773 containing only second 
and forth order terms. The second order term is given 

by 

f (% 2 +^)+/^2%, (45) 

where 

a = a+ (71 + y)^, = 0m- (46) 

The condition for a second order phase transition, that 
leads to finite values of 772 and 773 is a < \f3\. As we have 
a > and linear growth of \/3\ and quadratic growth 
of a — a with 771, the condition is fulfilled neither for 
large nor for small values of 771. But for intermediate 
values of 771 it can be fulfilled. This tendency back to- 
wards the original hexagonal symmetry in this or a sim- 
ilar form could be responsible for the appearance of ad- 
ditional Bragg peaks at the intermediate temperature of 
80-100K in Nao.sCoC^^ Note however, that it was spec- 
ulated that these Bragg peaks only exist over a narrow 
range of temperature. 



IX. CONCLUSIONS 

In this paper the properties of a high-symmetry multi- 
orbital model for the Co02-layer in combination with lo- 
cal Coulomb interaction are discussed. The tight-binding 
model is a zeroth order approximation to the kinetic en- 
ergy, as it only includes the most relevant hopping pro- 
cesses using Co-0 7r-hybridization. Nevertheless it pro- 
duces the hole pocket with predominantly a\ Q character 
around the T point, in agreement with both LDA calcu- 
lations and ARPES experiments. Furthermore, the three 
further pockets around the M points, although not seen 
in ARPES experiments, suggest that additional degrees 
of freedom that can not be captured in a single-band pic- 
ture could be relevant. The existence of identical hole, 
pockets in the BZ does however not produce pronounced 
nesting features. 



The local Coulomb repulsion of the ^-orbitals can be 
taken into account by an effective interaction of fermions 
with four different flavors, associated with the four hole 
pockets or the four inter-penetrating Kagome lattices. 
This effective interaction has a large discrete symmetry 
group, that allows to classify the spin- and charge-density 
operators, and to the determine for every mode the cor- 
responding coupling constant. 

It turns out that with an effective trigonal distortion, 
that splits the degeneracy between the T and the M 
points, general corrections to the quadratic part of the 
Hamiltonian, such as trigonal distortion or additional 
hopping terms, can be taken into account, provided they 
are small. This effective trigonal distortion reduces the 
symmetry of the Hamiltonian down to the space-group 
symmetries of the Co02-plane, by breaking the gauge 
symmetries of the effective interaction. 

Most coupling constants are negative for reasonable as- 
sumptions on the Coulomb integrals U, U' ', J' and Jh, 
but the ferromagnetic coupling constant is most negative 
and constitutes the dominant correlation. The charge 
and spin density wave instabilities without trigonal dis- 
tortion are easily described in a mean-field picture. In 
reciprocal space the degenerate bands split, and if bands 
belonging to different pockets hybridize, the BZ is folded. 
In real space different types of orderings are possible. The 
occupancy of the different t^g orbitals on different sites 
can be nonuniform, resulting in a charge ordering with 
nonuniform charge distribution on the Co-sites. Further, 
certain real or complex linear combinations of ti g orbitals 
can be preferably occupied on certain sites. In this case, 
the charge is uniformly distributed on the sites, but de- 
pending on the linear combinations of the orbitals, cer- 
tain space group symmetries are broken. The complex 
linear combinations of t<i g orbitals have in general a non- 
vanishing expectation value of the orbital angular mo- 
mentum. 

The tendency to these rather exotic states turns out 
to be smaller than the ferromagnetic tendency, and this 
dominance of the ferromagnetic state is even more en- 
hanced by the trigonal distortion. This is in good agree- 
ment with experiments, where ferromagnetic in-plane 
fluctuations have been observed by neutron scattering 
measurements in Nao.75Co02^ 3,24 There are also several 
reports of a phase transition in Nao.75Co02 at 22K to a 
static magnetic order, which is probably ferromagnetic 
in-plane but antiferromagnetic along the axis ^22^ 

In Nao.5Co02, a periodic Na-superstructure couples 
directly to a charge pattern in our model and crystallizes 
already at room temperature, whereas simple x \/3- 
superstructures, that would correspond x = 1/3 or x = 
2/3 do not couple. 

For general values of x the disordered Na-ions provide 
a random potential that couples to the charge degrees of 
freedom. Due to the incommensurability, this does not 
lead to long range order, but the short range correlations 
will also be influenced by the charge degrees of freedom 
in the C0O2 layers. This interaction between the Na- 
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correlations and the charge degrees of freedom could be 
the origin of the charge ordering phenomena at room 
temperature and the observation of inequivalent Co-sites 
in NMR experiments 

The overall agreement of our model with the experi- 
mental situation is good. Ferromagnetic fluctuations are 
dominant in our model and in experiments. Further- 
more, our model is based on a metallic state and allows 
for charge ordering and spin density ordering transitions 
without changing the metallic character of the state. Fi- 
nally, the clear Na superstructures, that were found at 
x = 0.5, can be understood quite naturally in this model. 

On the other hand there are still many open questions 
for the cobalt ates. Mainly the origin and the symmetry 
of the superconducting state of Na^CoC^-yHbO is still 
under debate. Unfortunately the Na content x = 0.35 
is beyond the validity of the approximations made in 
the derivation of our model. But also the samples with 
x > 0.5 have still many intriguing properties like the 
strongly anisotropic magnetic susceptibility, which shows 
the unusual Curie- Weiss temperature dependence. A 
possible description of the anisotropy of the magnetic 
susceptibility could be achieved by introducing a spin- 
orbit term into the kinetic energy. 

We hope that our model will be useful for further 
understanding of the rich experimental situation of the 
cobaltates. 
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APPENDIX B 



In this appendix, we provide some details concerning 
the derivation of the effective Hamiltonian in Eq. (|19|) . 
It is convenient to treat each term in Eq. (|T5|) separately. 
Let us start with the Hund's coupling. 
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The sum over the momenta in (|B2j) is restricted such that 
k + k 7 — q — q' equals a reciprocal lattice vector. (|B3|) 
follows from (|B2|) by using the definition of the pocket 
operators in Eq. (|BJ). The sum over the pocket indices is 
again restricted such that + Bj -f b& + B/ equals a recip- 
rocal lattice vector, whereas the sum over the momenta in 
the reduced BZ is simplified to an unrestricted sum over 
three momenta. Note that this simplification is valid for 
small pockets, because all the processes at the Fermi en- 
ergy are kept. (Small pockets means here 4k f < |bi|, 
this corresponds to a doping with x > 0.55) The next 
step is to go from orbital operators to the band operators. 
Restricting ourselves to the top band and taking into ac- 
count Eq. ([THjl we can simply substitute bli m(T — > -^bH a . 
Now we can sum over the orbital indices in Eq. (|B3|) and 
taking into account that the sum over the pocket indices 
is restricted we obtain the sum 
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APPENDIX A 



The equivalence of the two definitions for the "pocket- 
operators" made in Eq. (0) and in Eq. ((SJ) follows from 
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and for the Hund's coupling term 
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The restriction of the sum can be dropped, if we replace 
(ASik-1) with (26 ijk i-el jkl -6ii6 j k-Si j 6ki-\-36ikS j i). The 
terms proportional to Jh in the interaction of Eq. (|19[) are 
now obtained by dividing Eq. (|B5|) into two equal parts, 
rewrite one directly in terms of density density opera- 
tors, and rewrite the other in terms of density-density 
and spin-density spin-density operators using the SU (2) 
relation 25 a s5p 7 = 6 ' ai 6 : ps -+cr ai -cr -ps- Terms which renor- 
malize the chemical potential are dropped. All the other 
terms in Eq. (jTSjl are treated in the same way. 



The diagonal form of the tight-binding Hamiltonian in 
Eq. (|TT|) follows directly from the relation 
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mm — *Bj-(a m — a m /) mm 
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(A2) 



The symmetry group G of iJ e ff is a finite subgroup 
of U (4) that is generated by t the permutation matrices 
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V G S4 and the diagonal orthogonal matrices V G (Z2) 4 . 
G is a semi-direct product of 64 and the normal subgroup 
(Z2) 4 , this allows us to find the irreducible representa- 
tions of G, cf.51 The elements can be written in a unique 
way as (V,V) with V G <S 4 and V G (Z 2 ) 4 . The product 
of two elements (V,V)o (V ', V) is given by (V oV' V"). 
From this follows that if {V^V) is conjugate to (V 1 \T)'\ 

V is conjugate to V and the class of (V,V) G G can 
be labelled by the class of V G <S 4 . The elements of 
64 can be classified by writing them as disjunct cyclic 
permutations. We label the five classes as follows: e=l, 
f=(ab), g=(ab)(cd), h=(abc), i=(abcd). In total there 
are twenty 20 classes in G. The character table is shown 
in TABLE The character corresponding to the nat- 
ural representation of G by orthogonal 4x4 matrices is 
Xn- The representation on the 16 dimensional space V 
spanned by Q 0-15 , that was defined in section IIVI acts 
irreducibly on the subspaces V°, V 1-3 , V 4-9 and y 10-15 



with the characters xo? X7? X15 an d Xi6? respectively. 

With help of Schur's Lemma, it is now easy to show 
that the interaction H e ^ in the Basis Q 0-15 is diagonal, 
i.e. 

Qji^-fklQkl = Srr'K' S - (CI) 

and that the coupling constant A^ s depend only on the 
irreducible subspace. 

As discussed in section|3 the subgroup (Z 2 ) A describes 
gauge-symmetries, that are broken in the real system 
whereas the subgroup S4 describes the space-group sym- 
metries. The subgroup £4 consists of the classes ei, /1, 
gi, hi and i\. The irreducible representations of G are in 
general reducible for the subgroup S4. For example we 
have X7 = T 5 , X15 = r x © T 3 © T 5 and xie = T 4 © T 5 . 
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TABLE V: The character table for the symmetry group G of the effective Hamiltonian iJ e ff- The first line labels the classes 
and gives the number of elements in each class. The letters of the classes indicate classes of the subgroup <S 4 : e=l, f—{ab)^ 
g=(ab)(cd)j h=(abc)j i=(abcd). The characters appearing in our effective Hamiltonian are xi for Q°, X7 for Q 1-3 , xi5 for Q 4-9 
(real matrices) and xi6 for Q 10-15 (imaginary matrices), xn is the natural representation of G defined in section llVl The 
last column gives the reduction of the representations into irreducible representations of the subgroup <S 4 , that consists of the 
classes ei, /i, gi, hi and i\. 



